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Abstract  —  The  processing  of  remotely  measured  color 
imagery  involves  knowledge  of  the  diffuse  reflectance  of  the 
oceanic  water.  This  paper  presents  two  approaches  to  this 
problem.  The  first  one  is  analytical  and  consists  of  an 
approximate  solution  to  the  radiative  transfer  equation.  As  a 
result  it  gives  an  equation  for  diffuse  reflectance  of  the  sea  as 
a  function  of  inherent  optical  properties,  the  sun  elevation 
angle  and  the  ratio  of  the  direct  illumination  by  sun  to  the 
diffuse  illumination  by  sky.  The  second  approach  is  a 
numerical  Monte  Carlo  simulation.  The  results  of  this 
simulation  are  processed  to  produce  regressions  that  connect 
diffuse  reflectance  of  the  seawater  for  different  hydrooptical 
situations  to  the  sun  elevation  angle. 

INTRODUCTION 

In  papers  by  Haltrin  [1],  Gordon  [2],  Morel  and  Gentili 
[3]  it  was  shown  that  the  dependence  of  the  diffuse  reflection 
on  the  sun  zenith  angle  is  important  and  should  be 
incorporated  into  remote  sensing  algorithms. 

Two  approaches  to  this  problem  are  considered.  The  first 
approach  is  an  analytical  one  and  consists  of  approximate 
solution  to  the  radiative  transfer  equation.  As  a  result  it  gives 
Eqn.  (20)  for  diffuse  reflectance  of  the  sea  as  a  function  of  the 
inherent  optical  properties,  sun  elevation  angle  and  the  ratio 
of  the  direct  illumination  by  the  sun  to  the  diffuse 
illumination  by  the  sky.  The  second  approach  is  a  numerical 
Monte  Carlo  simulation.  It  is  based  on  the  modified  by  the 
author  [4]  approach  by  Kirk  [5].  The  results  of  this 
simulation  are  processed  to  produce  a  number  of  regressions 
that  connect  the  diffuse  reflectance  of  the  seawater  for  different 
hydrooptical  situations  with  the  sun  elevation  angle.  The 
results  of  both  methods  show  that  the  corrections  to  the 
diffuse  reflection  of  seawater  due  to  the  conditions  of 
illumination  can  reach  40%  at  certain  angular  elevations  and 
inherent  optical  properties  of  seawater. 

THEORETICAL  APPROACH 

Let  us  consider  a  homogeneous  sea  illuminated  by  the 
skylight  and  the  light  of  the  sun  elevated  at  degrees. 
According  to  the  Snellius  law,  the  direct  sunlight  enters  the 
sea  at  the  angle  cos"^  from  the  vertical  axis  Oz ,  directed 
from  the  sea  surface  to  the  sea  bottom,  where 

,  (1) 

here  is  a  refraction  coefficient  of  seawater. 

Let  us  start  with  the  system  of  two  flow  equations  for 
the  downward  Ej  and  upward  irradiances  proposed  in  Ref 
[1]  and  incorporated  later  into  the  elastic  part  of  Ref  [6]: 


‘  +  (2- fJ,){a-\-b^) 


Idz 

-(2-Ji)b,E,(z)^ 


E/z)  “  (2  +  ^) (z)  =  /(z) , 


L  dz 


(2) 


Here  z  is  the  depth  coordinate,  Jl  is  the  average  cosine  over 
the  irradiance  angular  distribution  in  the  sea  depth: 

,  EcOf. 

^  =  — ^^(3) 


1-coQ-hBcoQ 


I  l-l-2.x  +  -^.r(4-f’  5x)  '  <3  +  ^5 

b^  =  bB  is  the  backscattering  coefficient,  b  is  the  scattering 
coefficient,  R  =  is  the  probability  of 

backscattering,  p(j})  is  the  scattering  phase  function,  7}  is 
the  scattering  angle,  x  is  the  Gordon’s  parameter,  (O^-bi c 
is  the  single  scattering  albedo,  c  =  (3  +  /7  is  the  attenuation 
coefficient,  a  is  the  absorption  coefficient.  The  source 
functions  in  the  right  parts  of  Eqns.  (2)  are  equal  to: 

f{,z)  =  b^Est\^[-azl  B,),  a  =  a  +  2bg,  (4) 
where  ji^E^  is  the  sun  irradiance  just  below  the  sea  surface. 
In  Eqns.  (2)  the  direct  sunlight  is  taken  into  account  in  the 
form  of  the  source  functions  /(z). 


SOLUTIONS  FOR  IRRADIANCES 


Solutions  of  Eqns.  (2)  with  /  =  0  for  the  case  of  purely 
diffuse  illumination  of  the  sea  surface  are  given  in  Refs.  [1]. 
Let  us  find  the  solutions  of  these  equations  for  the  case  of 
combined  illumination,  i.  e.  for  the  case  when  /  ^  0.  Let  us 
consider  the  irradiance  of  the  skylight  penetrated  into  the  sea 
Eq  by  the  following  boundary  condition: 

=  (5) 

In  addition  let  us  accept  that  just  below  the  sea  surface  the 
directed  irradiance  of  the  sun  q  times  stronger  than  the 
irradiance  of  the  sky,  L  e.,  E^-qE^.  Let  us  also  define  the 
following  quantities:  let  E^{z)  be  the  portion  of  the  diffuse 
light  that  is  originated  from  the  scattering  of  the  skylight 
penetrated  into  the  ocean,  and  let  El  (z)  be  the  portion  of  the 
diffuse  light  that  is  originated  from  the  scattering  of  the 
sunlight  penetrated  into  the  sea. 

Let  us  introduce  the  following  definitions: 

/?„  =  E^f(0)/£,(0)  (6) 

is  the  diffuse  reflectance  of  the  infinitely  optically  deep  ocean 
illuminated  by  diffuse  light;  let 

R,  =  El{,0)l{^l^E,)  (7) 
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be  the  diffuse  reflectance  of  the  infinitely  optically  deep  ocean 
illuminated  by  directed  light  of  the  sun;  let 

+  (8) 

be  the  diffuse  reflectance  of  the  infinitely  optically  deep  ocean 
illuminated  by  the  combined  light  of  sun  and  sky.  Then  the 
diffuse  reflectance  of  the  sea  illuminated  by  the  natural  light 


’■  R^,  Sf^  - 


\  +  lu^qR,/R^ 


Let  us  look  for  a  solution  of  the  system  of  Eqns.  (2)  in  the 
form  of  the  sum  of  general  and  partial  solutions: 

here  -a^^-alp  is  a  negative  eigenvalue  of  Eqns  (2). 
Inserting  Eqns.  (10)  into  Eqns.  (2)  and  applying  boundary 
condition  (5)  we  have  following  equations  for  irradiances: 
Ej(z^  =  E°iz)  +  E^,iz),  £„(z)  =  £f(z)  +  £„"(z),  (11) 

E^iz)  =  E,e-^-\  E^{z)  =  KE,e-^-\  (12) 


EjU)  =  M^o 


Rq{^R^+Rs] 
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;  ^  Ep  Rf)  R, 
l-E,R^ 


a-jir 


^,=Jl-cos^hJnl-,  (19) 
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and  3)  the  diffuse  reflection  coefficient  of  the  sea  illuminated 
by  the  combined  light  of  the  sky  and  the  sun: 

^  _  (20) 

The  diffuse  reflection  coefficient  R^  depends  on  the  inherent 
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DIFFUSE  REFLECTANCE  COEFFICIENTS 

After  some  algebra  we  have  the  following  equations  for  the 
diffuse  reflection  coefficients  of  the  sea  illuminated  by  the 
different  kinds  of  light:  1)  the  diffuse  reflection  coefficient  of 
the  sea  illuminated  by  the  light  of  the  sky: 

Ji=  I  I.  (18) 

[1-juJ  ^]  l  +  2x  +  ^Jx(4  ^i-  5x) 

where  jc  =  b^/ia  +  bg)  =  B(D^I{\-o)^  +  BcOq)', 

2)  the  diffuse  reflection  coefficient  of  the  sea  illuminated  by 
the  direct  sunlight: 


Fig.  1.  Two-dimensional  density  plot  of  the  ratio  S^  as  a 
function  of  the  Gordon  parameter  x  and  the  Sun 
elevation  angle. 

optical  properties  of  the  water  b^  and  a,  the  sun  elevation 
angle  and  the  parameter  q  that  determines  the  ratio  of  the 
direct  sunlight  to  the  light  of  the  sky.  Correct  evaluation  of 
this  parameter  should  involve  solution  of  the  radiative  transfer 
problem  in  atmosphere  [7]  and  estimation  of  the  transmission 
through  the  air-sea  interface.  More  simple  approach  was 
proposed  by  Jerlov  [8].  According  to  Ref.  [8],  in  the  case  of 
the  clear  sky  the  ratio  q  depends  only  on  the  sun  elevation 
angle  h^.  Our  approximation  of  the  data  published  in  Ref.  [8] 
gives  the  following  formula: 

(7  =  0.25(I-h0.3/2j,  (21) 

here  is  measured  in  degrees.  Taking  into  account  Eqns. 
(18)-(21),  we  have  the  following  expression  for  the  ratio  of 
the  diffuse  reflectances  =  R^l  Rj. 

1  -  (1  -I-  -f  77(4  -  )  .^1  -  cos^  h,  /  nl 

S„  =  1 - IjA .  . . J .  (22) 

1  -I-  4/ (1  -f  0.3  /? J  y  1  -  cos^  /?,  /  nl 
Figure  1  shows  a  two-dimensional  density  plot  of  the  ratio 
Sf^  as  a  function  of  the  Gordon  parameter  x  and  the  sun 
elevation  angle  h^.  The  difference  between  maximum  and 
minimum  values  of  this  ratio  lies  in  the  range  of  40%. 

MONTE  CARLO  APPROACH 

This  approach  is  explained  in  detail  in  Ref.  [4].  The  main 
result  of  this  approach  related  to  the  influence  of  sun 
elevation  is  the  following  regression  equation: 

=  R,(9Q°)/RM, )  =  -^0  +  "  ^2  R-l  (23) 
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that  relates  the  ratio  of  diffuse  reflectance  coefficients  to  the 
cosine  of  the  sun  penetration  angle  .  The  coefficients  of 
these  regressions  are  given  in  Table  1.  The  difference  between 
the  theoretical  diffuse  reflectance  (20)  and  the  numerically 
simulated  values  do  not  exceed  20%. 


Table  1.  Coefficients  c^yC^,  and  C2  for  each  of  the  fifteen 
Petzold  phase  functions  [9,  10]. 


# 

^0 

Cl 

^2 

01 

0.03444 

0.288392 

1.58743 

0.29775 

0.999 

02 

0.01408 

0.565779 

3.23327 

1.66143 

0.995 

03 

0.01322 

0.576934 

3.18952 

1.60944 

0.997 

04 

0.01936 

0.967853 

3.86569 

1.89692 

0.999 

05 

0.01566 

1.024818 

4.03269 

2.00333 

0.999 

06 

0.08194 

0.852371 

3.32589 

1.46969 

0.999 

07 

0.09064 

0.810208 

3.24149 

1.42623 

0.999 

08 

0.14786 

0.477916 

2.38229 

0.90046 

0.999 

09 

0.01139 

-0.02428 

2.01621 

1.03560 

0.987 

10 

0.06660 

0.847081 

3.28663 

1.43518 

0.999 

11 

0.01265 

1.108715 

4.21634 

2.10450 

0.998 

12 

0.08947 

0.762611 

3.03790 

1.26988 

0.999 

13 

0.03550 

1.057443 

3.92840 

1.86609 

0.999 

14 

0.02131 

0.718718 

3.41994 

1.69663 

0.999 

15 

0.01355 

0.040812 

2.19756 

1.14971 

0.990 

CONCLUSIONS 

This  paper  presents  two  approaches  to  the  problem  of 
angular  dependence  of  the  sea  diffuse  reflection  coefficient. 
One  of  the  approaches  is  analytical  and  it  consists  of  the 
approximate  solution  to  the  radiative  transfer  equation.  The 
result  of  this  approach  is  Eqn.  (20)  for  diffuse  reflectance  of 
the  sea  as  a  function  of  inherent  optical  properties,  the  sun 
elevation  angle  and  the  ratio  of  direct  illumination  by  the  sun 
to  the  diffuse  illumination  by  the  sky.  The  second  approach  is 
a  numerical  Monte  Carlo  simulation.  The  results  of  this 
simulation  are  regressions  given  by  Eqn.  (23)  that  connect 
ratios  of  the  diffuse  reflectances  of  the  seawater  at  90°  and 
arbitrary  degrees  with  the  sun  elevation  angle. 

The  results  of  both  methods  show  that  the  corrections  to 
the  diffuse  reflection  of  seawater  due  to  the  illumination 
conditions  can  reach  40%.  This  effect  is  significant  and 
should  be  taken  into  account  in  processing  optical  remote 
sensing  data. 
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